Interference of a thermal Tonks gas on a ring 
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A nonzero temperature generalization of the Fermi-Bose mapping theorem is used to study the 
exact quantum statistical dynamics of a one-dimensional gas of impenetrable bosons on a ring. We 
investigate the interference produced when an initially trapped gas localized on one side of the ring 
is released, split via an optical-dipole grating, and recombined on the other side of the ring. Nonzero 
temperature is shown not to be a limitation to obtaining high visibility fringes. 
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A fundamental assumption at the heart of current 
proposals to realize integrated atom sensors is that the 
guided atom wavepackets will display interference phe- 
nomena when they are split, propagated, and recom- 
bined. On the other hand it is well known that the 
coherence properties of atomic gases are affected by di- 
mensionality, there being no true off-diagonal-long-range 
order in less than three dimensions, and this raises the 
issue of interference in restricted geometries. Motivated 
by recent theoretical arguments |l], |^| demonstrating that 
several stimulated processes for matter waves such as 
four-wave mixing, superradiance, and matter-wave am- 
plification can be achieved in degenerate fermion gases as 
well as in Bose-condensed gases, we have recently shown 
H H that a one-dimensional gas of hard-core bosons, or 
Tonks gas, at zero temperature can exhibit high visibil- 
ity interference fringes. These results suggest that the 
same mechanisms might be capable of overcoming the 
weakening of interference due to thermal excitation. The 
quantum Tonks gas is realized [|[ |(3) in a regime essen- 
tially opposite from that required for BEC, namely, the 
regime of low temperatures and densities and large pos- 
itive scattering lengths where the transverse mode be- 
comes frozen and the many-body Schrodinger dynamics 
becomes exactly soluble via a Fermi-Bose mapping the- 
orem 0, g, g, fu^ fy], |12). We shall extend that the- 
orem so as to obtain exact results for a model atom 
interferometer using a thermal Tonks gas. In particu- 
lar, we apply our results to the interference produced 
when an initially trapped gas localized on one side of 
a ring is released, split via an optical-dipole grating, 
and recombined on the other side of the ring. Such a 
study is currently of relevance due to experimental ef- 
forts |£| HI, HI lH [jj, H H, |tj| to fabricate atomic 
waveguides for matter wave interferometers. 

Model: The model consists of a ID gas of No hard 
core bosonic atoms on a ring. This situation can be re- 
alized physically using a toroidal trap of high aspect ra- 
tio R = 2L/£q where 2L is the torus circumference and 
£q the transverse oscillator length £q = ^/h/mujo with 
loo the frequency of transverse oscillations, assumed to 
be harmonic. The longitudinal (circumferential) motion 
can be described by a ID coordinate x along the ring 



with periodic boundary conditions. The pure-state quan- 
tum dynamics of the system is described by the time- 
dependent many-body Schrodinger equation (TDMBSE) 
ihd^B/dt = H^S> b with Hamiltonian 



H 



2m 



3 = 1 



dxj 



,x N ;t). 
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Here Xj is the ID position of the j particle, 
V'-B^ii ' " i^N'i t) is the N-particle wave function with pe- 
riodic boundary conditions, 



±2L---x N ;t) = 



■x N :t), 



(2) 

which is also symmetric under exchange of any two par- 
ticle coordinates, as is the many-body potential V. We 
consider the case of impenetrable two-particle interac- 
tions, the so-called Tonks-gas regime |, |, and this 
is conveniently treated as a constraint on allowed wave 
functions: 



if 



Xk 



1 < j < k < N, (3) 



rather than as an infinite contribution to V, which then 
consists of all other (finite) interactions and external po- 
tentials. 

We assume preparation of the system such that its 
initial state at time t = is expressed by a statistical 
density operator p = J2 a w a\^Ba(0))(^Ba(0)\ where 
the \^Ba(0)} are a complete set of orthonormal many- 
boson states with label a that stands for a set of quantum 
numbers, and the w a are nonnegative statistical weights 
summing to unity. Then the statistical average of any 
observable O at any later time t > is 



(4) 



where the states \^Ba(t)) evolve according to the 
TDMBSE with Hamiltonian (JTJ) . 

Statistical Fermi-Bose mapping theorem: The Fermi- 
Bose mapping theorem for pure quantum states ]7], ||, [[], 
[lC], |fl], is a mapping |* B ) = A\^ F ) from the Hilbcrt 
space of one-dimensional many-fermion states \^f) to 



2 



that of one-dimensional many-boson states \^b), hold- 
ing if the interparticle interaction contains a hard core 
and the Hamiltonian is of the form (Q). In Schrodinger 
representation the mapping operator A on TV-particle 
states consists of multiplication by the "unit antisym- 
metric function" 



A(xi,- ■ ■ ,x N ) 



Yl sgn(> fe 

l<j<k<N 



Xj), 



(5) 



where sgn(x) is the algebraic sign of x 
Xj, i.e., it is +!(-!) if x > 0(x < 0). 



Xk - 

Since 

A is constant except at nodes of the wavefunctions, 
the mapping converts antisymmetric fermionic solu- 
tions $>f(xi, ■ ■ ■ , %N', t) of the TDMBSE into symmetric 
bosonic solutions ^b{xi, ■ ■ ■ , xn', t) satisfying the same 
TDMBSE with the same Hamiltonian |l]), and satisfy- 
ing the same boundary conditions and constraint (||). In 
the Olshanii limit j5| (low density, tightly confining atom 
waveguide, large scattering length) the dynamics reduces 
to that of the impenetrable point Bose gas, which is 
then satisfied trivially (typ vanishing when any Xj — xe) 
due to antisymmetry. If the potential in Eq. (P is a 
sum of one-body external potentials, V(x±, ■ ■ ■ , x^; t) — 
Y^i* V(xi, t), the solutions of the fermion TDMBSE can 
be written as Slater determinants 0, |j, ^ ^ [□], 



1 



N 

det (f> ni (xj,t) 



(6) 



where rij G {0, 1, • ■ •} and the <p ni are orthonormal solu- 
tions of the single particle time-dependent Schrodinger 
equation (TDSE) for potential V(x, t). 

The generalization to quantum statistical evolution (^) 
of many-boson systems is now straightforward: 

(0(t)) B = ^w a {9 Pa (t)\A- l dA\9 Fa (t)) (7) 



where A -1 is the inverse mapping from the Hilbert space 
of bosonic states to that of fermionic states. This is 
most useful if the observable O commutes with the map- 
ping operator A in which case the identity A^ 1 A = 1 
implies that (0(t)) B — (0(t))p, reducing the statisti- 
cal evolution of the O in the impenetrable point Bose 
gas to that of the same observable in an ensemble, with 
the same statistical weights, of ideal Fermi gas states 
I^Fa) related to the corresponding impenetrable point 
Bose gas states $b ) by the pure-state mapping theo- 
rem [§11,0 

Single particle density: Interference fringes arising 
from spatial density modulation are due to first-order co- 
herence manifested in the single-particle density p(x,t); 
the corresponding operator commutes with A since it de- 
pends only on particle positions. It follows that the den- 
sity profiles for the Fermi and Bose problems are iden- 
tical, as can also be seen directly from the fact that 
A 2 (xi, ■ ■ ■ ,x N ) = 1 and hence |* B | 2 = |*f| 2 - 



The pure state mapping theorem ensures that corre- 
sponding Bose and Fermi states have the same parti- 
cle number N and energy levels E V (N), so that both 
systems have the same chemical potential /i, and there- 
fore corresponding states labeled by a = {N, E V (N)} 
will have the same grand canonical statistical weights 

W a = e i^-E u {N)) g|, Thug wc describe the gag 

of impenetrable bosons by a grand canonical ensemble 
of ideal Fermi gas states typ a (xx, ■ ■ ■ , £jv;£), assuming 
that the mean number of atoms No is sufficiently large 
to justify a grand canonical picture. At time t — we 
consider a basis of energy eigenstates in the external po- 
tential V(x, 0). The thermal average of the single particle 
density of the system is 

(o( x t)) = ^N=ol^uPN,EAx,t)e 

\f J ( x ' 1 ')/ ^oo e (nN-E v (N)) ' 

where the denominator is the grand partition function 
and pNE v (x,t) — 53i=i \4>ni(x,t)\ 2 is the spatial density 
of a TV-particle energy eigenstate of a free Fermi gas , its 
energy E V {N) = J2i=i e rn being the sum of the energies 
e ni of the single particle states <p ni (x,t). Some manipu- 
lation of the sums leads to the familiar form 



(9) 



where /„ = l/(e( e ™~^)/ fcf3T + 1) are Fermi-Dirac occu- 
pation numbers. The chemical potential fi is determined 
by SS^Lo fn — N . The density at a later time t > is 
obtained from Eq. (^) after propagating each orbital <j> n 
separately according to the TDSE for potential V(x,t). 

Initial condition: We assume that the ring has been 
loaded (say, by optical tweezers p!q]) with N atoms 
at temperature T, and for t < they are confined to 
a narrow segment of the ring by a trapping potential 
V(x, t < 0) assumed harmonic, with natural frequency 
to, over the spatial extent of the initial trapped gas. The 
grand canonical ensemble in Eqs. (H) and (§) describing 
the state of the atoms at t = is therefore characterized 
by the single particle energies e n = (n+ ^)hiu. The basic 
configuration is shown schematically in Fig. |f|(a). 

In order to discuss the time-evolution of the system, we 
designate the normal modes of a class of ID harmonic os- 
cillators by u n (x— x, w) where x are their mean positions 
and the parameter w relates their widths xgy/l + w 2 to 
the width of the initial trap xq — Wh/rnuj. Thus the 
modes of the initial trap potential are u n {x— x, w = 0). 
We choose the circumferential coordinates to have the 
trap center at x — L = —L. On unwrapping the ring 
about x — in Fig. |l|(b), the initial Hermite-Gaussian 
orbitals are split into two parts at the ends of the funda- 
mental periodicity cell [—L,L], and can be written as 

<f> n (x,t = 0) =u n (x + L,0) + u n (x- L,0) . (10) 

The asymptotic dependence of the harmonic oscilla- 
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FIG. 1: (a) Initial density of N hard core bosons trapped 
on a ring of circumference 2L. (b) By unfolding the ring we 
describe the system using a ID coordinate x £ [—L,L]. The 
coincident point x = — L = L is chosen at the center of the 
initial trapped gas. 



tor modes u n (x) ~ cos(\/2nx / xq) suggests the defini- 
tion of a critical wavevector for a thermal distribution, 
k = \/2^/xq, which would give a measure of the char- 
acteristic rate of expansion of atoms released from the 
initial longitudinal trap. We take it to correspond to 
fn = /o/2, so that for temperatures large compared to 
the Fermi temperature (Tp) the critical wavevector coin- 
cides with the thermal wavevector k rs yj2ksT /huj / x$, 
and for T — > we get the wavevector corresponding to 
the highest occupied mode k — > y / 2No/xQ. If the initially 
trapped gas is allowed to expand freely along the ring, 
the critical wave-vector gives a measure of the time to 
wrap around the ring t wrap = L/(hk/m). 

Some necessary conditions need to be satisfied in order 
that the initial atom cloud be accurately represented by 
a Tonks gas j|, |). First, the longitudinal energy must be 
small compared with the transverse excitation energy: at 
zero temperature this requires Nftw <C Hujq and at finite 
temperatures fcgT Hujq. Second, for the initial gas to 
be accurately described as an impenetrable gas of bosons 
we require fc|airj| <C 1, i.e. ksT <C muj 2 a 2 /2 ||. 

Optical- dipole grating: In order to produce interference 
from the initial trapped gas we turn off the harmonic trap 
at t = and apply a temporally short but intense spa- 
tially periodic potential of wavevector k. This spatially 
periodic grating may be produced over the spatial extent 
of the trapped gas, for example, using intersecting and 
off-resonant pulsed laser beams to produce an optical- 
dipole grating whose wavevector may be tuned by vary- 
ing the intersection angle. The applied periodic poten- 
tial then produces counter-propagating scattered atomic 
waves, or daughter waves, from the initial gas, or mother, 
with momenta ±fifc and these recombine on the opposite 
side of the ring at a time t r = L/(hk/m) = (k/k)t wrap . 
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FIG. 2: Scaled density profiles xop(x,t r ) for iVo = 100 atoms, 
n = 0.5, L/xq — 100 grating wavevector xok = 30 shown at 
the time of recombination t r for temperatures (a) T = 5 nK, 
(b) 20 nK and (c) 45 nK. The insets show the details of the 
fringes around the region of detection x = 0. 

infinite line. 

Following Rojo et. al. [|| we use a delta-function ap- 
proximation for the short pulse excitation of the periodic 
grating at t = 0, and implement the optical-dipole grat- 
ing using a phase-imprinting scheme according to which 
each orbital just after the pulse at t = + is changed to 



Ma:,0 + ) 



^ir/ cos[k(x- 
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T 







ik(x-L) 



-ik(x-L) 



(ii) 

(x,0), 



where rj is the grating amplitude. The second line as- 
sumes rj < 1 in which case the optical-dipole grating pre- 
dominantly produces two scattered copies of the initial 
mode travelling to the left and right with wavevectors ±fc 
in addition to the initial parent mode <f> n (x, 0). 

Interference on a ring: We wish to examine the pres- 
ence or absence of interference at the opposite side of the 
ring at time t = t r and how this depends on the temper- 
ature T and the grating vector k. After the grating is ap- 
plied the trapping potential is turned off, V(x, t > 0) = 0, 
so the subsequent quantum dynamics of the system is 
traced by freely propagating each orbital <j) n {x, + ). Us- 
ing the time evolved orbitals <p n (x,t r ) in Eq. (^|) we ob- 
tain the atom density at time t T as a sum of three terms 
p(x,t r ) = p m (x,t r ) + Pd(x,t r ) + p md (x,t r ), with 



For times * < 2t wrap the mother packet has not yet en- pm ^^ ^ = V /„ [u 2 n (x+L, ujt r ) + u 2 n {x-L : 



circled the ring and the gas expands freely as if on an 
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FIG. 3: Estimate of the visibility in the region of detection, as 
a function of temperature for three values of the optical-dipole 
grating wavevector (a) xok — 20, (b) 30 and (c) 45. 



+2u n (x+L, Lot r )u n {x—L 1 u>t r ) cos 



2xLu>t r 



pd{x,t r ) 



OC 



-[1 + cos(2kx)] } f n u n (x,tot r ). 



71=0 



(12) 



The first term p m is due to the expanding mother packet 
overlapping with itself, p d is due to the overlap of the 
daughters, and p m( i arises from the overlap of the mother 
and the daughter packets. In general, for t wrap > t r , 
which is desirable so that the mother packet does not 
wrap at the same time as the daughters recombine, p m d is 
small compared to the other terms: This term is present 
in our simulation below but we have refrained from re- 
producing the expression for it above. 

As a concrete example we consider No = 100 sodium 
atoms on a ring of circumference 2L — 100xo released 
from an initial harmonic trap of frequency lj — 2n x 4 
Hz. This corresponds to an oscillator length of x = 10.4 
pm and Fermi temperature Tp — 19.1 nK. In general, 
we find that for temperatures and applied wavevectors 
such that k <C k minimal interference fringes result when 
the daughter packets recombine at time t r . Recalling 
that the wrap time for the mother packet is given by 
twrap — (k/k)t r , then for k <C k the mother packet wraps 
around significantly by the time the daughters recombine, 
so p m provides a pedestal for interference fringes due to 
the overlap of the daughters pd, leading to reduced fringe 
visibility ||. In contrast for k >k pronounced interfer- 
ence fringes can appear when the daughter packets are 
recombined. This is illustrated in Fig. ^ where we fix the 
wavevector k of the optical grating and vary the temper- 
ature or alternatively k. Unit visibility fringes are seen 
in Fig. ||(a) where T = 5 nK, whereas reduced visibility 
fringes appear in Figs. |(b,c) where T = 20,45 nK. 

By noting that around x = the fringe pattern max- 
imum of p(x,t r ) is ~ p m (0,t r ) + pd(0,t r ), and the min- 
imum ~ p m (0,t r ), we obtain a measure of the fringe 
visibility as 



r = 



Pd(0,t r 



2p m (0,t r ) + p d (0,t r ) 



(13) 



This fringe visibility is plotted in Fig (||) as a function of 
temperature for three values of the grating vector k, with 
the clear and physically obvious trend that higher tem- 
peratures require a higher grating wavevector to obtain 



high visibility interference fringes. Physically there are of 
course limitations: For sodium atoms, the condition for 
a Tonks-gas isT^; (2x 10~ 5 x u ) 2 nK, for a transverse 
trapping frequency luq = 2'kvq Hz, and current experi- 
mental limits are about vq ~ 2 x 10 4 Hz p3] |. For the 
sodium yellow lines k ~ 10 m _1 , so for our parameters 
xok ~ 100, in keeping with the parameters used above. 

In conclusion, we have shown that even at nonzero 
temperature a strongly interacting ID gas of impenetra- 
ble bosons can show high visibility interference fringes on 
a ring. This remarkable result is of importance for cur- 
rent schemes to realize integrated atom interferometers 
in that it shows that neither many-body interactions nor 
nonzero temperature are fundamental limitations, in ID 
at least. Our results followed from a generalization of 
the zero temperature Fermi-Bose mapping which maps 
the strongly interacting ID boson problem to a ID gas 
of free fermions, and it is an open problem whether our 
conclusions can be extended to 2D and 3D interferome- 
ters. 
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